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The stability issue of Generalized modified gravitational models is discussed with par- 
ticular emphasis to de Sitter solutions. Two approaches are briefly presented. 
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1. The de Sitter stability issue 



It is well known that recent astrophysical data are in agreement with a universe 
in current phase of accelerated expansion, in contrast with the predictions of Ein- 
stein gravity in FRW space-time. Most part of energy contents, roughly 75 % in the 



^Jy universe is due to mysterious entity with negative pressure: Dark Energy. The sim- 

plest explanation is Einstein gravity plus a small positive cosmological constant. As 
an alternative, one may consider more drastic modification of General Relativity: 
Modified Gravity Models, see for exampleP^ 
^3 ' We shall deal with modified generalized models, described by a Lagragian density 

F(R,P,Q)® where R is the Ricci scalar, and P = R^R^ , and Q = R^pR^" 
are quadratic curvature invariants. In particular the Gauss-Bonnet topological in- 
variant reads G — R 2 — 4P + Q. The stability of the de Sitter solution, relevant for 
Dark energy, may be investigated in these Gauss-Bonnet models in several ways. 
We limit ourselves to the following two approaches: Perturbation of Esq. of Motion 
in the Jordan frame and Dynamical System Approach in FRW space-time. 



X 1 2. Stability of F(R, P, Q) model in the Jordan frame 

u : 

The starting point is the trace of the equations of motion, which is trivial in Einstein 
gravity R — —n 2 T, but, for a general F(R, P, Q) model, reads 

V 2 (3F R + RF' P ) + 2V fJ ,V v [(Fp + 2F' Q ) R^] - 2F + RF' R + 2(F' P + F' q ) = n 2 T . 

(1) 

Requiring R = Rq = Cts, Pq = Cts, and Qq = Cts one has de Sitter existence 
condition in vacuum 

2F - R F Ro - 2P F Po - 2Q' Qo = . (2) 

Perturbing around dS space, namely R = Ro + 5R, and with P = Pq + 5P, and 
Q = Qo + 5Q, observing that SP = ^-SR, and 5Q = ^-SR, one arrives at the 
perturbation Eq. 



- V 2 f5i? + M 2 5R = , 



(3) 
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in which the scalaron effective mass reads 



M 2 Ro Fk +^UlF' Po +F^) 

> \Rc[A R + A Q + A P + AHl (|F; o + F' Qo )} 



where 



Ar = F^ qRo + 6H*F>l Po + 4H 2 F£ oQa (5) 

A Q = 2Hl (F'> oPo + 6H*F£ oPo + 4if 2 F£ oQo ) (6) 

A P = ml (F^ Qo + 6H^ qPo + 4H^ oQa ) (7) 

Thus, if M 2 > 0, one has stability of the dS solution. In the particular case 
F(R, G), onehasPE] 

9F' 

Ro > 1 . (8) 



F 



In the case of a F(R) models, one has the well known condition Rg p°, > 1. 
3. Dynamical System Approach 

This approach has been used by many authors. One works in a cosmological setting, 
namely with a FRW metric, and the main idea consists in rewriting the general- 
ized Einstein-Friedman equations in an equivalent system of first order differential 
equations, introducing new dynamical variables fl; 

jti(t)=v(ti(t))- (9) 

Here the evolution parameter has been denoted by t. The critical (or fixed) points 
arc defined by v(Clo) — 0. The key point is: 

Hartman-Grobman theorem: The orbit structure of a dynamical system in 
the neighbourhood of a hyperbolic fixed point is topologically equivalent to the orbit 
structure of the associated linearized dynamical system, defined by a stability matrix 
M . 

Thus, in order to study the stability of the above non linear system of differential 
Eqs. at critical points, it is sufficient to investigate the related linear system of 
differential Eqs.: 

— 5fl(t) — MoSQ(t) , Mo Jacobian matrix evaluated at Qq (10) 

The solution of the linearization is well known and the evolution is determined by 
the signs of the eigenvalues of Mq. As a result, the non linear system is stable if all 
eigenvalues of the matrix Mq have negative real parts. 
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As an example, let us consider a modified model R + f(G). The related au- 
tonomous system in the two unknown quantities G and H reads 

The critical points are defined by G — and H = 0. Thus, we have the solutions 
24Hq = Go and Go/q — / o = 6Hq and these correspond to a de Sitter critical point 
with Gauss-Bonnet invariant. The linearized system around de Sitter critical point 
reads 

1 Xf 

SG = H SG-— T -SH, 5H=— m -AH 5H. (12) 
One can read off the stability matrix and the stability condition is „f fl , > 1, in 

J 

agreement with the previous approch. 



4. Concluding remarks 

Modified gravity may be seen as the phenomenological description of a fundamental 
unknown theory. From this point of view, corrections to Einstein-Hilbert action 
depending on higher order curvature invariants are likely to be expected (Lovelock 
gravity is an example). 

Among many existing approaches, two methods have been illustrated in order 
to investigate the stability of these models around de Sitter critical points, and the 
dS stability conditions has been derived within two possible approaches. 

These methods have owns advantages and problems, and, in our opinion, both 
permit to study critical points and stability for modified gravitational models de- 
pending on arbitrary geometric invariants, generalising the results obtained for F(R) 
models with other methods (see, for example^^) . We conclude noting that the dy- 
namical system approach has also been applied to non local F(R) mo dels Eili2J 
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